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Abstract 

Background: To obtain insights into Escherichia coli 0157:H7 (0157) survival mechanisms in the bovine rumen, we 
defined the growth characteristics and proteome of 0157 cultured in rumen fluid (RF; pH 6.0-7.2 and low volatile 
fatty acid content) obtained from rumen-fistulated cattle fed low protein content "maintenance diet" under diverse 
in vitro conditions. 

Results: Bottom-up proteomics (LC-MS/MS) of whole cell-lysates of 0157 cultured under anaerobic conditions in 
filter-sterilized RF (fRF; devoid of normal ruminal microbiota) and nutrient-depleted and filtered RF (dRF) resulted in 
an anaerobic 0157 fRF-and dRF-proteome comprising 35 proteins functionally associated with cell structure, motility, 
transport, metabolism and regulation, but interestingly, not with 0157 virulence. Shotgun proteomics-based analysis 
using isobaric tags for relative and absolute quantitation used to further study differential protein expression in 
unfiltered RF (uRF; RF containing normal rumen microbial flora) complemented these results. 

Conclusions: Our results indicate that in the rumen, the first anatomical compartment encountered by this human 
pathogen within the cattle gastrointestinal tract (GIT), 0157 initiates a program of specific gene expression that enables 
it to adapt to the in vivo environment, and successfully transit to its colonization sites in the bovine GIT. Further 
experiments in vitro using uRF from animals fed different diets and with additional 0157 strains, and in vivo using 
rumen-fistulated cattle will provide a comprehensive understanding of the adaptive mechanisms involved, and 
help direct evolution of novel modalities for blocking 0157 infection of cattle. 
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Background 

Escherichia coli 0157 (0157) have been implicated in 
several human outbreaks since their being established as 
foodborne pathogens in 1982; an estimated 63,153 illnesses, 
2,138 hospitalizations and 20 deaths occur annually in 
the United States [1-4]. Human disease ranges from 
self-limiting watery diarrhea to debilitating bloody diar- 
rhea that can advance into often fatal, extraintestinal, 
secondary sequelae in susceptible patients [3,4]. Cattle 
are the primary reservoirs for 0157, with their recto- 
anal junction (RAJ) serving as the colonization site at 
which these human foodborne pathogens persist [4,5]. 
However, the first bovine gastrointestinal tract (GIT) 
compartment that 0157 encounters is the rumen, where 
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the dynamic environment of regurgitating food particles, 
bio -fermentation, changing pH, and production of varying 
amounts of volatile fatty acids (VFA) poses challenges for 
0157 survival [6-8]. Studies thus far, have concentrated 
on the recovery of 0157 from the rumen, the in vitro 
0157 growth dynamics in modified rumen fluid or media 
with additives to mimic the rumen environment, expres- 
sion of select 0157 genes under controlled pH and VFA 
conditions, dietary effects on bacterial survival, and effects 
of select flora/metabolite on the growth/survival of 0157 
in the rumen or rumen fluid [6-11]. Despite this, however, 
a comprehensive study of the mechanisms used by 0157 
to survive the rumen environment is yet to be undertaken. 
Hence, as an initial step, we determined the repertoire 
of 0157 proteins (proteome) as expressed in vitro in 
harvested, rumen fluid (RF). We included RF of varying 
compositions (with and without normal flora, or depleted 
of nutrients essential for bacterial growth), with no ad- 
ditives, and used diverse culture conditions, to identify 
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bacterial factors that may enable 0157 adaptation to 
the rumen. 

Methods 

Bacterial strain, inoculum preparation and animals 

Wild-type 0157 strain 86-24 (Shiga toxin (Stx) 1 -negative, 
Stx 2-positive; motile; clinical isolate) was used in this 
study [12]. Overnight culture of 0157 in Luria-Bertani 
(LB) broth, grown at 39°C with aeration was used to pre- 
pare log-phase sub-cultures of the same in 50 ml LB broth, 
under the same growth conditions. Bacteria harvested 
from the log-phase cultures at an OD 600 0.5-0.6, washed 
and re-suspended in sterile 0.9% saline, were used to 
inoculate various rumen fluid (RF) or LB aliquots as 
described under 'Culture conditions and processing for 
proteomics'. All 0157 cultures were confirmed sero- 
logically using latex agglutination kits (Remel Inc., Lenexa, 
KS). Two rumen-fistulated Holstein cows, routinely used 
as rumen fluid 'donors' at the National Animal Disease 
Center (NADC, Ames, I A) with approval from the 
NAD C- Animal Care and Use Committee, were used in 
this study. Both animals, approximately 1 year of age, were 
fed the NADC Maintenance Diet (corn silage, grass hay, 
520 pellets, protein supplements) at 25% fiber and 10% 
protein, with ad-lib access to water through out. 

Unfiltered (uRF), Filtered (fRF), and Depleted RF (dRF) 

Rumen fluid samples collected from the two animals 
(Samples A and B; Tables 1 and 2), on separate days, 
were used to prepare the RF-preparations for each ex- 
periment set (Experiment I and II). Two liters of RF was 
collected 2-3 hr post-feeding to allow for rumination to 
occur, at each sampling time [10,13]. RF was strained 
through cheesecloth to remove large feed particles, and 
poured into collection flasks; pH was recorded on site 



Table 1 Biochemical characteristics of rumen fluid used 
to analyze growth patterns of 0157 strain 86-24 in 
Experiment I 



Sample 
analysis 


Depleted rumen fluid 


Filtered rumen fluid 


Sample A 


Sample B 


Sample A 


Sample B 


pH 1 


7.9 


7.6 


7.6 


7.7 




Volatile Fatty Acids (urn/ml; VFA) 




Total VFA 


324 


207 


211 


157 


Acetic acid 


201 


142 


144 


112 




(62%) 2 


(69%) 


(68%) 


(71%) 


Propionic acid 


41 


28 


31 


23 




(13%) 


(14%) 


(15%) 


(15%) 


Butyric acid 


43 


20 


16 


10 




(13%) 


(10%) 


(8%) 


(6%) 



Vh, post-depletion and/or post-filtration of the depleted and filtered rumen 
fluid samples, respectively. 

2 Percent individual volatile fatty acid of the total is shown in parenthesis. 



and an aliquot frozen at -80°C for volatile fatty acid 
(VFA) analysis. Approximately 500 ml of the strained RF 
was stored as the unfiltered RF (uRF) at 4°C. 

One half of the remaining strained RF was processed 
as follows to generate filtered RF (fRF). The strained RF 
was centrifuged at 27,000 x g for 30 mins at 18°C, at least 
3 times, to remove particulate matter and pressure fil- 
tered using a 0.5 \i pre-filter and a 0.2 \i filter in tandem 
(Pall Corporation, Port Washington, NY). The fRF was 
collected into sterile bottles and stored at 4°C after re- 
cording the pH and freezing an aliquot for VFA analysis. 

To prepare dRF, the other half of the remaining 
strained RF was first subjected to depletion, a process 
that involves exhaustion of residual nutrients in the RF 
by exploiting metabolic activities of the resident micro- 
flora, prior to the centrifugation-filtration steps. Specific- 
ally, the depletion process was initiated by adjusting the 
strained RF pH to 6.8-7.0, and incubating it under anaer- 
obic conditions, at 39°C for four days. The strained RF 
was held in flasks fitted with stoppers bearing valves to 
release the fermentation gases throughout the incuba- 
tion, following which the depleted RF was centrifuged 
and filtered as described above. This depletion protocol 
was adapted from previously described methods with no 
extraneous substrates added to the RF prior to depletion 
[11,14]. The pH of the resultant filter-sterilized dRF was 
recorded and aliquots set aside for VFA analysis prior to 
storage at 4°C in sterile bottles. 

pH and volatile fatty acids (VFA) analysis 

Initial rumen fluid pH measurements were taken during 
collection by using a portable pH meter (Thermo Fisher 
Scientific Inc., Waltham, MA) [8,11]. Subsequently, the 
pH meter or pH paper was used (pH range 5.0-8.0; 
Micro Essential Laboratory Inc., Brooklyn, NY), to record 
pH of the processed RF and media. VFA concentrations in 
rumen fluid and its preparations were determined by 
capillary gas chromatography of their butyl esters, 
as described previously [15,16], on an Agilent 6890 N 
gas chromatograph (Agilent Technologies, Inc., Santa 
Clara, CA). 

Culture conditions, and processing for proteomics 

RF preparations from Samples A and B were analyzed 
separately per experiment set, and each analysis in turn 
was conducted in duplicate. In Experiment I, 5 ml LB, 
dRF, or fRF media were aliquoted separately into 85, 
16 x 150 mm tubes. Of these, five tubes per media were 
used as uninoculated controls. The remaining 80 tubes 
were inoculated with 0157. To create anaerobic culture 
conditions, half of these tubes were transferred into the 
anaerobic Coy Chamber for 72 hrs, sealed and inocu- 
lated within the chamber and then removed. The log- 
phase 0157 culture, re-suspended in 0.9% saline was 
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Table 2 Biochemical characteristics of rumen fluid used to analyze growth patterns of 0157 strain 86-24 in 
Experiment II 

Sample Depleted rumen fluid Filtered rumen fluid Unfiltered rumen 

analysis Sample A Sample B Sample A Sample B Sample A 

pH 1 7.6 7.4 7.7 7.2 6.4 

Volatile Fatty Acids (um/ml; VFA) 



Total 


203 


205 


144 


153 


210 


Acetic acid 


139 


140 


103 


110 


141 




(68%) 2 


(68%) 


(72%) 


(72%) 


(67%) 


Propionic acid 


28 


28 


21 


23 


32 




(14%) 


(14%) 


(13%) 


(15%) 


(15%) 


Butyric acid 


19 


19 


9 


10 


20 




(9%) 


(9%) 


(6%) 


(7%) 


(10%) 



^H, post-depletion and/or post-filtration of the depleted and filtered rumen fluid 
2 Percent individual volatile fatty acid of the total is shown in parenthesis. 

inoculated to a starting OD 600 0.05-0.06, into all the 80 
tubes, which were then incubated at 39°C with shaking, 
along with the uninoculated control tubes. 0157 was 
grown to an OD 600 of 0.8-1.0, before harvesting cells 
from each tube by centrifugation at 7,000 rpm, 15 min 
at 4°C. Bacterial cells from like media, whether derived 
from RF-samples A or B, were pooled together and 
washed three times with an equal volume of ice-cold 
sterile phosphate buffered saline (PBS; pH 7.4), and 
processed to obtain cell lysate and pellet fractions for 
bottom-up proteomic analysis [17]. 

In Experiment II, uRF was included to the media (LB, 
dRF, fRF) being evaluated and aliquoted as described 
above. However, the 0157 inoculum diluted in saline to 
the starting OD 600 0.05-0.06 was placed in sterile dialysis 
tubing (Spectra/Por Type F, PVDF: 80,000 kDa cut off; 
Serva Electrophoresis, Heidelberg, Germany) and sus- 
pended within the uRF containing tubes [18]. This was 
to ease the recovery of 0 157 from the complex uRF mi- 
lieu and the colony counts recovered from the tubings 
matched those obtained by magnetic recovery of 0157 
from directly inoculated uRF (data not shown). 0157- 
innoculated LB, dRF, fRF, and uRF were incubated for 
48 h, anaerobically, before harvesting cells and process- 
ing for proteomic analysis [17] using iTRAQ. For this 
experiment, bacterial cells from like media were pooled 
together but kept separate between preparations derived 
from RF-samples A and B. The culture conditions used 
in Experiment II correlated with ruminal conditions and 
feed turnover rates [19-21]. In both experiments, OD 600 
of each tube was recorded relative to uninoculated control 
tubes, centrifuged at 10,000 rpm for 10 min to remove 
any sediments or particulate matter which could interfere 
with the spectrophotometer reading. In addition, pH, and 
colony counts (on LB agar) were determined from the five 
uninoculated and ten inoculated tubes at different time 
points, for comparison. 



i fluid 
Sample B 

6.7 

165 
104 
(63%) 

30 
(18%) 

17 
(10%) 
samples, respectively. 

Bottom-up proteomics using liquid chromatography 
tandem mass spectrometry (LC-MS/MS) 

This proteomic analysis was done at the Proteomics 
Division, ICBR, University of Florida, Gainesville, Florida. 
0157 cell pellet and lysate fractions from Experiment I 
(LB, dRF, fRF) were concentrated using spin filters (MW 
cutoff 5000 Daltons), and digested with trypsin prior to 
tandem mass spectrometry (MS/MS) as described pre- 
viously [17]. The enzymatically-digested samples were 
injected onto a capillary trap (LC Packings PepMap) 
and desalted for 5 min with a flow rate of 3 ul/min of 
0.1% v/v acetic acid. The samples were loaded onto an 
LC Packing® CI 8 Pep Map nanoflow HPLC column. 
The elution gradient of the HPLC column started at 3% 
solvent B, 97% solvent A and finished at 60% solvent B, 
40% solvent A for 95 min for protein identification. 

Solvent A consisted of 0.1% v/v acetic acid, 3% v/v 
acetonitrile (ACN), and 96.9% v/v H 2 0. Solvent B con- 
sisted of 0.1% v/v acetic acid, 96.9% v/v ACN, and 
3% v/v H 2 0. LC-MS/MS analysis was carried out on a 
hybrid quadrupole-TOF mass spectrometer (QSTAR elite, 
Applied Biosystems, Framingham, MA). The focusing 
potential and ion spray voltage was set to 225 V and 
2400 V, respectively. The information-dependent acquisi- 
tion (IDA) mode of operation was employed in which a 
survey scan from m/z 400-1800 was acquired followed by 
collision-induced dissociation (CID) of the four most in- 
tense ions. Survey and MS/MS spectra for each IDA cycle 
were accumulated for 1 and 3 s, respectively. 

Tandem mass spectra were extracted by ABI Analyst 
version 2.0. All MS/MS samples were analyzed using 
Mascot (Matrix Science, London, UK; version 2.2.2). 
Mascot was set up to search NCBI with taxonomy Bacteria 
database assuming the digestion enzyme trypsin. Mascot 
was searched with a fragment ion mass tolerance of 
0.50 Da and a parent ion tolerance of 0.50 Da. Iodoace- 
tamide derivative of Cys, deamidation of Asn and Gin, 



Kudva et al. BMC Microbiology 2014, 14:48 
http://www.biomedcentral.eom/1 471 -21 80/1 4/48 



Page 4 of 18 



oxidation of Met, were specified in Mascot as variable 
modifications. Scaffold (version Scaffold-03-3-2, Proteome 
Software Inc., Portland, OR) was used to validate MS/MS 
based peptide and protein identifications. Peptide identifi- 
cations were accepted if they could be established at 
greater than 95.0% probability as specified by the Peptide 
Prophet algorithm [22]. Protein identifications were ac- 
cepted if they could be established at greater than 99.0% 
probability and contained at least 2 identified unique pep- 
tides. Proteins with single peptide hits were included if 
they exhibited high confidence based on low false discov- 
ery rates [23]. Relative protein abundance was estimated 
using the normailized total spectral counts [24]. Protein 
probabilities were assigned using the Protein Prophet al- 
gorithm [25] . Proteins that contained similar peptides and 
could not be differentiated based on MS/MS analysis 
alone were grouped to satisfy the principles of parsimony. 

Quantitative proteomics using isobaric tags for relative 
and absolute quantification (iTRAQ) 

0157 cell pellet and lysate fractions from Experiment II 
(LB, dRF, fRF, uRF; cultured 48 h, anaerobically) were 
analyzed by iTRAQ. Samples were processed, trypsin 
digested, and labeled with various iTRAQ reagents as 
described earlier [26], in accordance with the manufac- 
tures instructions for the iTRAQ 4-plex kit (Amine- 
Modifying Labeling Reagents for Multiplexed Relative 
and Absolute Protein Quantitation, Applied Biosystems, 
Foster City CA). Labeled peptides were combined, dried 
in one tube, and held at -80°C until use. A modification 
of the previously used protocol was used to analyze 
these labeled peptides that were resuspended in mobile 
phase A (72 mM triethlyamine in H20, pH 10 with 
acetic acid) at a concentration of 200 [ig/[il and incu- 
bated for 1 hour in a sonic-water bath at RT. 100 ug of 
sample was injected into a Waters 1525 u Binary HPLC 
(Waters Corporation, Milford, MA) with a Waters XBridge 
CI 8, 3.5um, 1 x 100 mm column in mobile phase A and 
ran isocratically for 6 minutes. The gradient consisted of, 
0-20% mobile phase B (72 mM triethlyamine in ACN, 
52 mM acetic acid), over 34 minutes; 20-40% over 
20 minutes; and finally 40-100% over 2 minutes, at a 
flow rate of 100 ul/ minute throughout the entire gradi- 
ent [27]. Two-minute fractions were collected, dried in 
a vacuum centrifuge, and resuspended in nano-HPLC 
buffer A (95% H 2 0: 5% ACN and 0.1% formic acid). Based 
on previous experience we combined, 3 fractions before 
and after, the fractions that contained the majority of the 
eluted peptides. 

Fractions from the first dimension chromatography 
were injected on a second dimension of chromatography 
using a Proxeon Easy-nLC (Thermo Fisher Scientific, 
West Palm Beach, FL) connected to the mass spectrometer. 
The second dimension chromatography used a trapping 



column (Proxeon Easy-Column, 2 cm, ID 100 um, 5um, 
120A, C18) and an analytical column (Proxeon Easy- 
Column, 10 cm, ID 75 um, 3 um, 120A, C18). The gradi- 
ent using a mobile phase A (95% H20: 5% acetonitrile 
and 0.1% formic acid) and mobile phase B (5% H20: 95% 
acetonitrile and 0.1% formic acid). The gradient was, 0% B 
for 3 minutes, 0%-8% B from 3-5 minutes, 8-18% B 
from 5-85 minutes, 18-30% B from 85-100 minutes, 
30-90% B from 100-105 minutes, and held at 90% B 
from 105-120 minutes at continuous flow rate through- 
out the gradient of 300 nl/min. The analytical column 
was connected to a PicoTip Emitter (New Objectives, 
Woburn, MA; FS360-75-15-N-20) and together attached 
to a LTQ OrbiTrap Velos Pro (Thermo Fisher Scientific, 
West Palm Beach, FL) mass spectrometer using the 
Proxeon Nanospray Flex Ion Source. The capillary 
temperature was set at 275°C and spray voltage was 
2.9 kV. The mass spectrometer was used in a data 
dependent method. In MS mode, the instrument was 
set to scan 300-2000 m/z with a resolution of 30,000 
FWHM. A minimal signal of 20,000 could trigger tan- 
dem MS and 10 consecutive MS/MS were possible. 
High-energy collision-induced dissociation (HCD) was 
used to resolve the iTRAQ reporter ions, 113-117. 
The normalized collision energy was set to 35 and repeat 
mass exclusion was set to 120 seconds. 

Tandem mass spectra were extracted and charge state 
deconvoluted by Proteome Discoverer version 1.4. Charge 
state deconvolution and deisotoping was not performed. 
All MS/MS samples were analyzed using Mascot, Sequest 
(XCorr Only; Thermo Fisher Scientific, San Jose, CA, 
USA; version 1.3.0.339) and X! Tandem (GPM.org; version 
CYCLONE (2010.12.01.1)) assuming digestion with tryp- 
sin. A custom E. coli database was generated by combin- 
ing the fasta files from uniprot.org from the following 
E. coli strains: 12009/EHEC, 2009EL-2050, 2009EL- 
2071, 2011C-3493, 11128/EHEC, 0157:H7, EC4115/ 
EHEC, TW14359/EHEC, and 11368/EHEC. This E. coli 
fasta file consists of 47,819 entries and was generated 
in May 2013. Mascot, Sequest (XCorr Only) and X! 
Tandem were searched with a fragment ion mass toler- 
ance of 0.100 Da and a parent ion tolerance of 10.0 
PPM; carbamidomethyl of cysteine and iTRAQ4plex of 
lysine and the n-terminus were specified as fixed modi- 
fications while deamidation of asparagine and glutamine, 
oxidation of methionine and iTRAQ4plex of tyrosine 
were specified as variable modifications. Scaffold (version 
Scaffold_4-0.6) wa s used to validate MS/MS based pep- 
tide and protein identifications, as described above for 
'Bottom-up Proteomics'. The 0157-proteome as expressed 
in LB was used as the reference against which all the other 
0157-proteomes were compared. Two biological replicate 
samples (Sample A and B), corresponding to the duplicate 
experiments described under 'Culture conditions, and 
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processing for proteomics' above, were analyzed separ- 
ately. In addition, each sample was analyzed twice (Run 
A and Run B; technical replicates) to cover the entire 
spectra of proteins in these samples. Only proteins that 
were consistently identified were selected for analysis. 

Statistics and bioinformatics 

The Student t-Test (two-tailed) was used to evaluate dif- 
ferences between the means of the 0157 optical densities 
and viable counts recovered from the different cultures 
and a values of p < 0.05 was considered significant. Puta- 
tive functions were determined by querying the Conserved 
Domain Database (CDD) at http://www.ncbi.nlm.nih.gov/ 
Structure/cdd/wrpsb.cgi, and associated metabolic path- 
ways were determined using the KEGG pathway database 
at http://www.genome.jp/kegg/pathway.html. Cellular and 
sub-cellular locations of proteins were determined as de- 
scribed previously [17]. 

Results 

pH and VFA content 

The pH and VFA concentrations were comparable 
amongst all rumen fluid samples, indicating consistency 
in maintenance diet being fed and the ruminal chemistry 
between the two animals enrolled in the study (Tables 1 
and 2). The pH of the uRF ranged from 6.4-6.7 at collec- 
tion [28-31] but attained a more neutral pH after filtering, 
as seen with dRF (pH 7.4-7.9) and fRF (pH, 7.2-7.7) in 
both experiments (Tables 1 and 2). Concentrations of 
three VFAs, acetate, propionate and butyrate, were closely 
analyzed as these vary the most with changes in the forage 
versus starch compositions of the feed, and are of rele- 
vance to both host and bacterial growth. Consistent with 
the 25% forage and 10% protein diet that these cattle were 
being fed, the RF comprised a higher percentage of acetate 
[28-31]. Acetate ranged from 72-62%, compared to the 
13-18% propionate and 6-13% butyrate concentrations 
across the uRF, dRF and fRF samples in both experiments, 
irrespective of procedures used to prepare dRF and fRF 
(Tables 1 and 2). LB broth (pH 7.0-7.2) did not contain 
added VFAs. 

0157 growth characteristics 

Log phase 0157 cultures, set up for the two experi- 
ments, were at 0.5-0.6 OD 600 , respectively, with viable 
counts around 1 x 10 8 cfu/ml. Hence, when each medium 
was inoculated to a starting 0.05-0.06 OD 600 , the corre- 
sponding 0157 counts were at -1-5 x 10 7 cfu/ml. In both 
experiments, 0157 grew to an OD 600 of 1.0 within 2 h in 
LB media, aerobically and anaerobically as anticipated, 
with an increase in viable count to 4 x 10 8 cfu/ml and 
the final culture pH at 6.0-6.2. However, significant dif- 
ferences were observed between aerobic and anaerobic 



growth patterns of 0 157 when cultured in dRF, fRF and 
uRF preparations. 

In Experiment I, 0157 cultured in dRF and fRF 
achieved an average OD 600 of 0.6-1.0 in 48 h aerobically, 
but remained at a low OD 600 of <0.2 anaerobically, even 
after 14 days of incubation. Irrespective of the ODs, 
viable 0157 was recovered from all cultures, but the 
viable counts at 10 6 (dRF)-2 x 10 7 (fRF) cfu/ml aerobic- 
ally, and at 10 5 (dRF)-2 x 10 5 (fRF) cfu/ml anaerobically 
(data not shown) appeared to be static or decreasing. The 
pH for dRF and fRF cultures at the end of incubation was 
around 7.7 (aerobic)-7.3 (anaerobic). Similar 0157 growth 
results were observed upon anaerobic culture for 48 h in 
dRF, fRF and uRF, in Experiment II (Figure 1), with the pH 
for uRF cultures being 6.8 at end of incubation. This 
was despite these media being prepared with RF from a 
separate animal and a shorter anaerobic incubation 
period than in the first experiment, thereby verifying 
the observations made initially. Here, the cultures reached 
an average OD 600 of 0.97 (LB), -0.03 (dRF), -0.04 (fRF) 
and -0.03 (uRF) in 48 h, with 0157 viable counts of 
2 x 10 8 cfu/ml (LB), 4 x 10 5 cfu/ml (dRF), 3 x 10 6 cfu/ml 
(fRF) and 1 x 10 6 cfu/ml (uRF), respectively. 

Significant differences were observed among the op- 
tical densities and viable counts of LB cultures versus 
RF-preparation cultures, under all growth conditions. 
However, differences between the RF-preparations were not 
always significant (Figure 1). For instance, in Experiment II, 
the p values for the 0157 viable counts were: LB: dRF, 
p = 0.0379; LB: fRF, = 0.0385; LB: uRF, = 0.0381, 
dRF: fRF, p = 0.0121 and fRF: uRF, p = 0.0655; dRF: uRF, 
p = 0.1077. 

Proteomics analysis 
(i) Bottom-up 

LC-MS/MS analysis of the 0157 cell pellet and lysate 
fractions generated in Experiment I provided insights 
into the proteins being expressed by 0157 in different 
media, under different growth conditions and at ex- 
tended incubation time points. A total of 585 protein 
(2284 spectra) hits were identified by setting minimum 
characteristics for the identification confidence. How- 
ever, of these only 218 0157 proteins matched a higher 
threshold cut off, with 90% protein-80% peptide prob- 
ability in the Scaffold Viewer, and hence, were selected 
for analysis. The 218 0157 proteins were differentially 
expressed: 90 only under aerobic conditions, 37 only 
under anaerobic conditions and 91 under both condi- 
tions (data not shown), accounting for fewer proteins 
under anaerobic conditions. Interestingly, none of the 
0157 proteins expressed aerobically or anaerobically in 
either media were associated with direct virulence (e.g., 
the Locus of Enterocyte Effacement [LEE] -encoded pro- 
teins or Shiga toxins) but were primarily associated with 



Kudva et al. BMC Microbiology 2014, 14:48 
http://www.biomedcentral.eom/1 471 -21 80/1 4/48 



Page 6 of 18 



/XO.OOOI 




B 



p 0.0381 



4.0x10 s i 
3.0xl0«- 
2.0xlO«- 
1.0x10 s - 

4.0xl0 6 
3.0x10*- 
2.0x10* • 
1.0x10* 
6.0x10* 
4.0x10* 
2.0x10* 
0 



0 0385 




! 










ll 


iii 









Figure 1 Growth characteristics of 0157 in Experiment II, following anaerobic incubation for 48 h, in LB and RF-preparations. Optical 
densities (OD 600 ) and viable counts (colony forming units [cfu]/ml) # with the standard error of means, are shown in graph A and B, respectively. 
The p values shown on the graphs were calculated using the Student t-Test (significant, p < 0.05). 



sequences homologous to other E. coli genomes (Backbone) 
(Additional file 1: Table SI). Considering that the rumen 
is an anaerobic microbiome, the 128/218 0157 proteins 
expressed anaerobically were examined in greater de- 
tail These proteins were either unique to growth in LB 
(93/128), dRF (2/128), fRF (10/128) or, expressed in more 
than one media (14/128 in LB/dRF/fRF, 9/128 in dRF/fRF) 
(Figure 2). Specifically, the 35 proteins expressed anaer- 
obically in fRF and dRF (unique and shared combined), 
were functionally associated with the osmotic adaptation 
pathway (OsmE), anaerobic respiration and oxidative 
stress pathway (YggE, MoaB, DmsB, FdoH), heat stress 
response (HchA), carbon starvation response (Sip), en- 
ergy metabolism and biosynthetic pathways (glycolytic/ 
gluconeogenesis pathway, amino acid biosynthesis: AldoC, 
Crr, AnsB, PykF, Eno, GpmA, GadpH, CysK, Ttc, AhpC, 
YhcB), chaperones (DnaK, GroEL, HchA), transport 
(LamB, ManX, FadL, RbsB), outer membrane proteins/ 
porins/channel (OmpC, TolC, YdeN, Sip, OmpA), tellur- 
ite resistance (TerD), lysozyme inhibitor (Ivy), chemotaxis 
(GgbP), and motility (FliC) (Table 3; Additional file 1: 
Table SI). 

(ii) iTRAQ 

To more closely examine and quantify 0157 protein ex- 
pression in the bovine rumen, especially in the uRF, the 
anaerobic 0157-proteome expressed in LB, dRF, fRF and 
uRF after 48 h incubation was compared using iTRAQ, 
in Experiment II. Data generated in two runs for each 
biological replicate was condensed to create a single 
comprehensive file per sample, and the files for the two 
biological replicate samples compared (Additional file 2: 



Table S2) to identify unambiguous proteins. Using the 
anaerobic 0157-proteome expressed in LB as the refer- 
ence, a total of 394 0157 proteins that were either dif- 
ferentially or similarly expressed in dRF, fRF, and uRF 
were identified (Figure 3, Additional file 2: Table S2). Of 
the cumulative 35 0157 proteins expressed anaerobically 
in dRF and fRF, and identified via Bottom-up proteomics, 
10 were not identified using iTRAQ in the second experi- 
ment (Table 3). Overall, only 134 proteins were common 




Figure 2 Distribution of 128 anaerobically expressed 0157 

proteins, identified using bottom-up proteomics, amongst the 

media tested. LB, Luria-Bertani broth; dRF, depleted and filtered 

rumen fluid; fRF, filtered rumen fluid. 
v> . J 



Table 3 01 57-proteome expressed under anaerobic conditions in dRF and fRF in Experiment I 



Protein/Function/Pathway; Name Accession Molecular Number of Peptides 

Number Weight (kDa) 1 (Relative Abundance) 2 











dRF 


fRF 


2,3-bisphosphoglycerate-dependent phosphoglycerate mutase; GpmA 


gi 


157155502 


29 kDa 


1 (6) 


1 (4) 


3-isopropylmalate dehydrogenase/Amino acid Biosynthesis; TtC 3 


gi 


1 70684236 


40 kDa 


1 (6) 


1 (4) 


Alkyl hydroperoxide reductase protein C/Energy; AhpC 


gi 


15800320 


21 kDa 


0 


1 (4) 


Anaerobic dimethyl sulfoxide reductase/Anaerobic growth/Oxidative Stress; DmsB 


gi 


145756 


23 kDa 


1 (6) 


1 (4) 


Chain A, Crystal Structure Of The Long-Chain Fatty Acid Transporter^ransport; Fadl 


gi 


203282230 


47 kDa 


1 (6) 


1 (4) 


Chain A, Crystal Structure Of Ggbp. Glucose-galactose binding protein/Chemotaxis, Transport; Ggbp 3 


gi 


126030485 


33 kDa 


0 


1 (4) 


Chain A, Structure Of Ivy/ Lysozyme inhibitor; Ivy 


gi 


29726212 


15 kDa 


1 (6) 


1 (4) 


Chaperone protein, stabilizes proteins under heat stress/Heat Stress Related; HchA 


gi 


15802400 


31 kDa 


0 


1 (4) 


Chaperonin, type 1 protein/Protein folding^ransport; GroEL 3 


gi 


15834378 


57 kDa 


3(18) 


3(11) 


Cysteine synthase/Amino acid transport and Metabolism; CysK 


gi 


145686 


35 kDa 


1 (6) 


1 (4) 


Cytochrome d ubiquinol oxidase subunit Ill/Oxidative phosphorylation/Energy; YhCB 3 


gi 


157148804 


15 kDa 


0 


1 (4) 


D-ribose transporter subunit B^ransport; RbsB 


gi 


110644091 


31 kDa 


1 (6) 


0 


DNA-binding transcriptional activator/Osmotically-inducible lipoprotein E; OsmE 


gi 


15802150 


12 kDa 


2(18) 


2(7) 


DNA-directed RNA polymerase subunit alpha^ranscription; RpoA 


gi 


123444073 


37 kDa 


0 


1 (4) 


Flagellin/Flagellar assembly/Motility; FliC 


gi 


15802358 


60 kDa 


3(24) 


4(19) 


Formate dehydrogenase-O, iron-sulfur subunit, energy metabolism/Anaerobic Respiration, Glyoxylate & Dicarboxylate Metabolism; FdoH 


gi 


15804482 


33 kDa 


0 


1 (4) 


Fructose-bisphosphate aldolase/Glycolysis, Gluconeogenesis, Amino acid Biosynthesis; AldoC 


gi 


161984958 


38 kDa 


3(24) 


4(19) 


Glucose-specific PTS system component, phosphorylation^ransport; Crr 


gi 


15802950 


18 kDa 


0 


1 (4) 


Glyceraldehyde 3-Phosphate Dehydrogenase; GadpH 3 


gi 


1421424 


35 kDa 


2(12) 


1 (4) 


Hypothetical protein CKO_00658/ Uncharacterized; DedA 3 


gi 


157144929 


21 kDa 


0 


1 (4) 


Hypothetical protein EcE24377A_0553/Glyoxylate Utilization; GlxB 


gi 


157157046 


29 kDa 


1 (6) 


0 


Hypothetical protein ECP_291 1 /Oxidative Stress; YggE 


gi 


110643066 


25 kDa 


0 


2(7) 


L-asparaginase II, induced by anaerobiosis/Nitrogen and Amino acid Metabolism; AnsB 


gi 


157157301 


37 kDa 


0 


4(15) 


Maltoporin/Receptor for lambda phage^ransport; LamB 3 


gi 


110644375 


50 kDa 


1 (6) 


1 (4) 


Molecular chaperone/Protein folding^ransport; DnaK 


gi 


157159481 


69 kDa 


7(42) 


4(15) 


Molybdopterin biosynthesis protein B/Cofactor Biosynthesis/ Oxidative Stress; MoaB 


gi 


15800533 


19 kDa 


0 


1 (4) 


Outer membrane channel protein, efflux of hydrophobic molecules^ransport; TolC 


gi 


110643281 


54 kDa 


0 


2(7) 


Outer membrane porin protein C^ranport of small molecules/Osmotic; OmpC 3 


gi 


15802768 


41 kDa 


1 (6) 


5(22) 


Outer membrane protein II, porin, receptor, integrity/Membrane Stability; OmpA 


gi 


146983 


26 kDa 


2(12) 


3(11) 


Outer membrane protein induced after carbon starvation, stationary phase, environmental stress/ Membrane stability; Sip 


gi 


110807343 


27 kDa 


0 


1 (4) 


Phosphopyruvate hydratase: enolase/Glycolysis, Gluconeogenesis; Eno 3 


gi 


15832893 


46 kDa 


0 


1 (4) 



Table 3 0157-proteome expressed under anaerobic conditions in dRF and fRF in Experiment I (Continued) 



PTS system, mannose-specific II AB component/phosphotransferase^ransport; ManX 


gi|1 10641934 


35 


kDa 


0 


1 (4) 


Putative sulfatase/lnorganic ion transport and metabolism^ran sport; YdeN 3 


gi|1 10641672 


63 


kDa 


0 


1 (4) 


Pyruvate kinase/Glycolysis, Gluconeogenesis, Amino acid Biosynthesis; PykF 


gi|1 10805653 


59 


kDa 


0 


1 (4) 


Tellurium resistance protein/Stress related; TerD 


gi|135596 


20 


kDa 


0 


1 (4) 



VDa, Kilodalton. 

2 Relative abundance based on normalized total spectral counts. 
3 Proteins not identified in Experiment II (see Table 4). 
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Figure 3 Log fold changes in the expression of 0157 proteins, identified using iTRAQ, in media tested under anaerobic conditions. The 

0157-proteome expressed in LB was the reference against which the regulation of 0157 proteins in other media was determined. The scatter 
plots represent 0157 proteins expressed in the context of the 155 up-regulated in uRF (Panel A), 82 down-regulated in uRF (Panel B) and 
157 with no change in expression levels in uRF (Panel C). LB, Luria-Bertani broth; dRF, depleted and filtered rumen fluid; fRF, filtered rumen 
fluid; uRF, unfiltered rumen fluid. 



to the results of the two experiments, indicative of 
incubation-time related differences in the number and 
type of proteins expressed. Differentially expressed 0157 
proteins in the iTRAQ dataset distributed as 298/394 in 
dRF (169, up-regulated, 129, down-regulated), 241/394 in 
fRF (162, up-regulated, 79, down- regulated) and 237/394 
in uRF (155, up-regulated, 82, down- regulated) (Table 4). 
Interestingly, similar expression patterns were observed 
between 0157 proteins expressed in dRF and uRF; 90% of 
dRF-differentially regulated and 71% dRF-no change pro- 
teins were similarly expressed in uRF. This may have been 
due to shared growth conditions (nutrient limitation)/ 
signals in these two media. The competing microflora 
in uRF may have decreased nutrients in that media. 

As observed with the Bottom-up proteomics results, 
none of the well-established 0157 virulence factors were 
identified in either media after 48 h of anaerobic incuba- 
tion (Additional file 2: Table S2). Specifically, the 155 
up-regulated (Table 4), uRF-0157 proteins could be 
functionally associated with osmotic adaptation (MdoG, 



CreC, OsmE, YjbJ), oxidative stress pathway (KatG, 
DmsA_ynfE, DmsA_ynfF, YggE), heat shock response 
(HdlD, GrpE), carbon starvation response (SspA, SspB), 
anaerobic respiration (HybC, ErpA, GlpE), pH adaptation/ 
acid resistance (SpeA, SpeB, Mrp), energy metabolism: 
degradation of carbon compounds (GalT, BglX, EbgA, 
MtlD), glycolysis/ gluconeogenesis (GpmA, SucB, FdhO, 
FbaB, GloB, NuoC, AceF, PoxB, SdhA), amino acid me- 
tabolism (GcvT, GcvP, HslV, IlvE, GlnA, TnaA), nitrogen 
and glycerolipid metabolism (GlpK), DNA degradation 
(RecB), biosynthetic pathways: fatty acid (FabG, AccD, 
AccB), amino acids (DapB, ArgH, AsD, IlvA), nucleotides 
(PyrI, PyrC, PurH, GlmU, CmK, GuaA, GuaB, GmK,), cel- 
lulose (BcsC), cofactors/carriers (MoaC, GshB), vitamins 
(RibB, PdxJ), chaperones (fimbrial usher protein,), trans- 
port (HmuV, FadL, PlaP, MacB, OppA, NikA, SecD, 
ManZ, PotA, YajC, EtpE), storage (BfR, FtnA), multi- 
drug efflux systems (AcrA, MdtC), tellurite resistance 
(TerC), serine proteases (DegP, DegQ, EspP), outer mem- 
brane proteins/porins/channel (AsmA, LptE, Lpp, NagA, 



Table 4 0157 proteins up-regulated under anaerobic conditions in uRF, in Experiment II 






3" 7^ 

S O 


Protein/Function/Pathway; Name 


Accession Number 


Molecular Weight ( kDa) 1 


^ < 


2,3-bisphosphoglycerate-dependent phosphoglycerate mutase/Glycolysis/Energy; GpmA 


GPMA_EC057 


29 kDa 


% % 


3-oxoacyl-[acyl-carrier-protein] reductase/Fatty acid biosynthesis; FabG 


Q8X8I5_EC057 


26 kDa 


S S 
o ^ 


3,4-dihydroxy-2-butanone 4-phosphate synthase, Riboflavin (Vitamin B2) biosynthesis/ Cofactor Bisoynthesis; RibB 


RIBB_EC057 


23 kDa 


3 n 

fD 

Q_ ^ 


30S ribosomal protein S6/Protein Translation; RpsF 


RS6_EC057 


15 kDa 


icrobiology 2 
central.com/ 


4-hydroxy-tetrahydrodipicolinate reductase, Leucine synthesis/Ami no acid Biosynthesis; DapB 


DAPB_EC057 


29 kDa 


50S ribosomal protein L10/Protein Translation; RplJ 


RL10_ECO57 


18 kDa 


50S ribosomal protein L18/Protein Translation; RpIR 


RL18_EC057 


13 kDa 


014, 
147' 


Acetyl-coenzyme A carboxylase carboxyl transferase subunit beta/Fatty acid Biosynthesis; AccD 


ACCD_EC057 


33 kDa 




Acriflavine resistance protein A/Multidrug efflux system^ran sport; AcrA 


B5Z3X3_EC05E 


44 kDa 


00 £ 


ADP-L-glycero-D-manno-heptose-6-epimerase/LPS core Biosynthesis/Heat induced; HldD 


HLDD_EC057 


35 kDa 


00 


Agmatinase/Polyamine Biosynthesis/Acid Resistance; SpeB 


SPEB_EC057 (+1) 


34 kDa 


AidA-l adhesin-like protein/Adherence; AidA 


K0AT24_ECO1C 


141 kDa 




Aminomethyltransferase/Nitrogen, Amino Acid Metabolism; GcvT 


GCST_EC057 


40 kDa 




Argininosuccinate lyase, Arginine synthesis/ Amino Acid Biosynthesis; ArgH 


ARLY_EC057 


50 kDa 




AsmA protein/LPS Biogenesis; AsmA 


B5YUC4_EC05E (+5) 


69 kDa 




Aspartate carbamoyltransferase regulatory chain/Nucleotide Biosynthesis; Pyrl 


PYRLEC057 


17 kDa 




Aspartate-tRNA ligase/Aminoacyl-tRNA Biosynthesis; AspS 


SYD_EC057 (+1) 


66 kDa 




Aspartate-semialdehyde dehydrogenase/Amino acid Biosynthesis; AsD 


DHAS_EC057 


40 kDa 




ATP synthase subunit b/Energy Production; AtpF 


ATPF_EC057 


17 kDa 




ATP-dependent helicase/DNA Replication, Repair; HrpB 


Q8X904_ECO57 (+1) 


91 kDa 




ATP-dependent protease subunit/Proteolysis; HsIV 


HSLV_EC057 


19 kDa 




Bacterioferritin/lron storage and Transport; BfR 


Q8X890_ECO57 


18 kDa 




Bacteriophage replication gene A protein/Predicted Phage replication; ECH741 15_3058 


B5YUH3_EC05E 


86 kDa 




Beta-D-glucoside glucohydrolase, degradation of small carbon compunds/Biosyntheis of secondary metabolites; BglX 


Q8X668_EC057 


83 kDa 




Bifunctional N-acetylglucosamin-1-phosphate-uridyl transferase/Cell Wall Synthesis; GlmU 


GLMU_EC057 (+2) 


49 kDa 




Bifunctional purine biosynthesis protein/Purine Biosynthesis; PurH 


PUR9_EC057 


57 kDa 




Biofilm regulator/Biofilms, Adherence; BssR 


BSSR_EC057 


15 kDa 




Biosynthetic arginine decarboxylase/Polyamine Biosynthesis, Acid Resistance; SpeA 


SPEA_EC057 


74 kDa 




Biotin carboxyl carrier protein of acetyl-CoA carboxylase/Fatty Acid Biosynthesis; AccB 


BCCP_EC057 (+8) 


17 kDa 




Branched-chain-amino-acid aminotransferase/Amino acid Transport and Metabolism; llvE 


ILVE_EC057 


34 kDa 




Catabolite repression sensor kinase for PhoB alternative sensor for pho regulon/Global Regulaor; CreC 


Q8XB07_ECO57 


52 kDa 


to 
a> 


Catalase-peroxidase 1 /Prevent Cell, DNA damage/Oxidative Stress; KatGI 


KATG1_EC057 


80 kDa 


o 
o 

— h 


Cell Division protein; ZapB 


ZAPB_EC05E (+8) 


9 kDa 


00 



Table 4 0157 proteins up-regulated under anaerobic conditions in uRF, in Experiment II (Continued) 

Cellulose synthase subunit/Biofilms, Adherence; BcsC 

Chromosome partition protein/Cell Division; MukB 

Curli production assembly/transport subunit/Biofilms, Adherence; CsgG 

Cyclic pyranopterin monophosphate synthase accessory protein/Cofactor Biosynthesis; MoaC 

Cytidylate kinase/Nucleotide Biosynthesis; Cmk 

Dihyd rolipoyllysine-residue succinyltransferase component of 2-oxoglutarate dehydrogenase 

Complex/Energy Metabolism; SucB 

Dihydroorotase/Nucleotide Biosynthesis; PyrC 

DNA helicase/DNA Replication, Transcription; RecB 

DNA invertase from prophage CP-933H/DNA Replication; PinH 

DNA topoisomerase IV subunit B/DNA Replication; ParE 

DNA-directed RNA polymerase subunit alpha/DNA Transcription; RpoA 

DNA-directed RNA polymerase subunit beta/DNA Transcription; RpoB 

DNA-directed RNA polymerase subunit omega/DNA Transcription; RpoZ 

Elongation factor P-like protein^ranslation, Protein synthesis; YeiP 

Evolved beta-D-galactosidase alpha subunit/Degradation of small carbon compounds; EbgA 
Ferritin-1/lron Uptake and Storage; FtnA 

Fimbrial usher family protein/Chaperone^ransport; ECH741 15_21 18 

Formate dehydrogenase-0 major subunit/Energy Metabolism, Anaerobic Respiration; FdoH 

Fructose-bisphosphate aldolase, class 1 /Glycolysis, Gluconeogensis; FbaB 

Galactose-1 -phosphate uridylyltransferase/Galactose Metabolism; GalT 

Glucans biosynthesis protein G/Osmotic Adaptation; MdoG 

Glutamate-tRNA ligase/Amino acyl tRNA synthesis; GltX 

Glutamine synthetase/Amino acid Biosynthesis; GlnA 

Glutathione synthetase/Cofactors, carriers Biosynthesis; GshB 

Glycerol kinase/Glycerolipid Metabolism; GIpK 

Glycine dehydrogenase [decarboxylating]/Amino acid Metabolism; GcvP 
GMP synthase [glutamine-hydrolyzing]/Nucleotide Biosynthesis; GuaA 
GTPase-activating protein^ranscriptional Activator; Yihl 
Guanylate kinase/Nucleotide Biosynthesis; GmK 
Hemin import ATP-binding protein^ran sport; HmuV 
Histidine-tRNA ligase/Amino acyl tRNA synthesis; HisS 

HTH-type transcriptional regulator/Maltooligosaccharide Uptake and Metabolism; MalT 



MUKB_EC057 (+2) 170 kDa fa 

B5YVQ8_EC05E (+1) 31 kDa o | 

3 o 

MOAC_EC057 1 7 kDa 2_ ^ 

n n' 

KCY_EC057 25 kDa £ | 

{D_ 6' 

ft 

0D02 EC057 44 kDa < nj 

— — » o 

PYRC_EC057 39 kDa § £ 

Q8X6M9_EC057 134 kDa £g 

o oo 

Q8X7L1_EC057 22 kDa 

Q8XBQ4_EC057 70 kDa 

RPOA_EC057 37 kDa 

RPOB_EC057 151 kDa 

RPOZ_EC057 10 kDa 

EFPL_EC057 22 kDa 

Q8XAM9_EC057 119 kDa 

FTNA_EC057 19 kDa 

B5Z1W3_EC05E (+5) 81 kDa 

Q7A9A6_EC057 (+1) 113 kDa 

B5YV44_EC05E (+1) 38 kDa 

Q8X943_EC057 40 kDa 

OPGG_EC05E 59 kDa 

SYE_EC057 54 kDa 

GLNA_EC057 52 kDa 

GSHB_EC057 35 kDa 

GLPK_EC057 56 kDa 

GCSP_EC057 104 kDa 

GUAA_EC057 59 kDa 

YIHLEC057 19 kDa 

KGUA_EC057 (+2) 24 kDa 

HMUV_EC057 29 kDa 

SYH_EC057 47 kDa 

MALT_EC057 103 kDa 



Table 4 0157 proteins up-regulated under anaerobic conditions in uRF, in Experiment II (Continued) 

Hydrogenase-2 large chain/Energy Metabolism/Anaerobic Respiration; HybC 

Hydroxyacylglutathione hydrolase/Pyruvate Metabolism; GloB 

lnosine-5'-monophosphate dehydrogenase/Nucleotide Biosynthesis; GuaB 

Iron-sulfur cluster insertion protein/Growth using alternate electron acceptors; ErpA 

Long-chain fatty acid transport protein^ran sport; FadL 

Low-affinity putrescine importer/Transport; PlaP 

LPS-assembly lipoprotein/LPS Biogenesis; LptE 

Macrolide export ATP-binding/permease protein^ransport; MacB 

Major outer membrane lipoprotein/Membrane Integrity; Lpp 

Mannitol-1 -phosphate 5-dehydrogenase/Fructose, Mannose Metabolism; MtlD 

Membrane, TerC family/CBS/transporter associated protein^ransport^ellurium Resistance; TerC 

Methyl-accepting chemotaxis protein III, ribose and galactose sensor receptor/Chemotaxis; Trg 

Minor curlin protein/Biofilms, Adherence; CsgB 

Multiple Resistance and pH adaptation protein, antiporter^ransport; Mrp 

Multidrug resistance efflux protein^ransport; MdtC 

Multiphosphoryl transfer protein/Fructose, Mannose Metabolism; FruB 

N-acetylglucosamine-6-phosphate deacetylase/LPS Biosynthesis; NagA 

NADH-quinone oxidoreductase subunit C/D/Electron transfer/Energy Metabolism; NuoC 

NifU-like protein/Energy Production, Conversion; NifU 

Non-LEE-encoded type III secreted effector; EspX7 

Oligopeptide ABC transporter, periplasmic oligopeptide-binding protein^ran sport; OppA 

Osmotically-inducible lipoprotein E/Global Regulator; OsmE 

Outer membrane lipoprotein/Cell Envelope Biogenesis; SlyB 

Outer membrane protein, porin, receptor, integrity/Membrane Stability; OmpA 

Outer membrane protein assembly factor/Membrane Assembly, Antibiotic Resistance; BamA/YaeT 

Outer membrane protein assembly factor/Membrane Assembly, Antibiotic Resistance; BamD/YfiO 

Outer membrane protein, efflux protein^ransport; TolC 

Outer membrane protein W/Outer Membrane Biogenesis; OmpW/YciD 

Oxidoreductase subunit/Oxidative Stress; DmsA_YnfF 

Penicillin-binding protein activator/Li pid binding/Putative Adhesin; LpoB/YcfM 

Peptidase B/Protein Metabolism; PepB 

Peptide chain release factor 1 /Protein Translation; PrfA 

Peptidyl-prolyl cis-trans isomerase A/Post-translational Modification; PpiA 



GL02_EC057 


28 kDa 


| 5. 


B5Z0X7_ECO5E 
ERPA_EC057 


54 kDa 
12 kDa 


BMC Microbiology 20' 
.biomedcentral.com/1 ' 


FADL_EC057 
PLAP_EC057 
LPTE_EC057 


49 kDa 

50 kDa 
21 kDa 


MACB_EC057 


71 kDa 




LPP_EC057 
MTLD_EC057 


8 kDa 
41 kDa 


14:48 

■2180/14/ 


B5YUC1_EC05E (+6) 


59 kDa 


00 


C8U899_ECO10 


59 kDa 




B5YVR3_EC05E (+1) 


17 kDa 




B5YV63_EC05E (+1) 


40 kDa 




C8TU07_ECO26 (+2) 


111 kDa 




PTFAH_EC057 


40 kDa 




NAGA_EC057 


41 kDa 




NUOCD_EC057 (+2) 


69 kDa 




NIFU_EC057 


14 kDa 




C6URC5_EC05T (+2) 


73 kDa 




B5YYE8_EC05E 


63 kDa 




OSME_EC057 


12 kDa 




SLYB_EC057 


16 kDa 




B5\T87_EC05E 


38 kDa 




BAMA_EC057 


91 kDa 




BAMD_EC057 


28 kDa 




B5YR81_EC05E (+2) 


54 kDa 




L)oXLDD_LLOb/ 


23 kDa 




C6UV29_EC05T (+6) 


90 kDa 




LPOB_EC057 


23 kDa 




PEPB_EC057 


46 kDa 
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RF1_EC057 


41 kDa 


PPIA_EC057 


20 kDa 



Table 4 0157 proteins up-regulated under anaerobic conditions in uRF, in Experiment II (Continued) 

Peptidyl-prolyl cis-trans isomerase/Post-translational Modification; FklB 
Periplasmic binding protein for nickel/Amino acid Transport and Metabolism; NikA 
Phenylalanine-tRNA ligase alpha subunit/Amino acyl tRNA synthesis; PheS 
Phenylalanine-tRNA ligase beta subunit/Amino acyl tRNA synthesis; PheT 
Poly (A) polymerase l/DNA Transcription; PcnB 
Proline-tRNA ligase/Amino acyl tRNA synthesis; ProS 
Protease IgAI, Serine protease/Protection; EspP 
Protein elaB/Uncharacterized; ElaB 

Protein grpE, prevents aggregation of denatured proteins/ Heat and Hyperosmotic Shock-Related; GrpE 
Protein translocase subunit^ransport; SecD 

PTS system, mannose-specific transporter subunit IID^ransport; ManZ 
Putative anaerobic dimethyl sulfoxide reductase chain A; DmsA_YnfE 

Putative DNA replication factor encoded within cryptic prophage CP-933P/Hypothetical; Z6069 
Putative endopeptidase of prophage CP-933X/Hypothetical; Z1877 
Putative carboxypeptidase/Hypothetical; YagX 

Putative lipoprotein induced during stationary phase/Stress Response; YbjP 

Putative homeobox protein/Regulator; YbgS 

Putative lipoprotein/Membrane protein; LppC 

Putative membrane protein, peptidase/Uncharacterized; YibP 

Putative multimodular enzyme/Energy Metabolism; Z3719 

Putative pectinesterase, localizes to cellular poles/Membrane protein; YbhC 

Putative replicase/DNA Replication; Z5187 

Pyridoxine 5'-phosphate synthase/Vitamin B6 (Pyridoxine) synthesis; PdxJ 

Pyruvate dehydrogenase (Dihydrolipoyltransacetylase component)/Energy Metabolism; AceF 

Pyruvate oxidase/Degradation of small carbon compunds; PoxB 

RNase E/RNA Degradation; RnE 

Serine protease/Protection; DegP 

Serine endoprotease/Protection; DegQ 

Single-stranded DNA-binding protein/DNA Replication; SsB 

Soluble cytochrome, electron transport/Energy Metabolism; CybC 

Spermidine/putrescine import ATP-binding protein^ransport; PotA 

Stringent starvation protein A, stationary phase induced acid tolerance/Global Regulator; SspA 
Stringent starvation protein B, CIpXP protease specificity enhancer/Global Regulator; SspB 



Q8X5U3_EC057 (+1) 


59 kDa 


| 5. 


SYFA_EC057 (+2) 
SYFB_EC057 


37 kDa 
87 kDa 


BMC Microbiology 20' 
.biomedcentral.com/1 ' 


PCNB_EC057 (+3) 

SYP_EC057 

K0AWD8_ECO1C 


54 kDa 
64 kDa 
146 kDa 


ELAB_EC057 


11 kDa 




C8U980_ECO10 (+2) 
SECD_EC057 


22 kDa 
67 kDa 


14:48 

■2180/14/ 


B5YQW0_ECO5E (+1) 


31 kDa 


00 


Q7ABM3_EC057 (+1) 


88 kDa 




Q8XAD9_EC057 


28 kDa 




Q8X704_ECO57 


12 kDa 




Q8X6I4_EC057 (+1) 


91 kDa 




Q8X6N7_EC057 (+1) 


19 kDa 




Q8X948_EC057 


13 kDa 




B5YUN1_EC05E 


20 kDa 




Q8XDE2_EC057 


47 kDa 




Q8XBF4_EC057 


82 kDa 




Q8X891_EC057 


46 kDa 




Q8XBZ7_EC057 


37 kDa 




PDXJ_EC057 


26 kDa 




Q8X966_EC057 


66 kDa 




Q8X6L4_EC057 


62 kDa 




Q8X8J5_EC057 


118 kDa 




B5Z0E1_ECO5E (+1) 


49 kDa 




QoXyh \_tLUb/ 


47 kDa 




SSB_EC057 


19 kDa 




C562_EC057 (+4) 


14 kDa 




POTA_EC057 


43 kDa 
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SSPA_EC057 


24 kDa 


SSPB_EC057 


18 kDa 



Table 4 0157 proteins up-regulated under anaerobic conditions in uRF, in Experiment II (Continued) 







Kud 
http 


DHSA_EC057 


64 kDa 


d < 
% % 


Q8X8J6_EC057 


30 kDa 




C8TJL4_EC026 (+2) 
Q8X467_EC057 


12 kDa 
56 kDa 


BMC Microbiology 20' 
.biomedcentral.com/1 ' 


NUSG_EC057 
CPXR_EC057 
C8TGU9_EC026 


21 kDa 
26 kDa 
72 kDa 


Q8XBF1_EC057 (+6) 


73 kDa 




IF2_EC057 


97 kDa 


14:48 

■2180/14/ 


TIG_EC057 (+1) 


48 kDa 


TNAA_EC057 


53 kDa 


00 


082884_EC057 


56 kDa 




WZC_EC057 


79 kDa 




Q8X3G9_EC057 (+1) 


17 kDa 




Q8X7H8_EC057 


14 kDa 




Q8X2Z1_EC057 


10 kDa 




YGGE_EC057 


27 kDa 




YQJD_EC057 


11 kDa 




YHBJ_EC057 


32 kDa 




YAJC_EC057 


12 kDa 




YJBJ_EC057 


8 kDa 




YEJL_EC057 


8 kDa 




PYRH_EC057 


26 kDa 




SYV_EC057 


108 kDa 





Succinate dehydrogenase flavoprotein subunit/Energy Metabolism; SdhA 

Tat-linked quality control protein/DNAse activity; TatD 

Thiosulfate sulfurtransferase/Anaerobic Respiration, Energy Metabolism; GIpE 

Threonine deaminase (Dehydratase)/Amino acid Biosynthesis; MvA 

Transcription termination/antitermination protein/Modulates DNA Transcription; NusG 

Transcriptional regulatory protein/Envelope Stress Response Protein/Downregulates LEE; CpxR 

Transketolase 1, thiamin-binding/Non-oxidative Metabloism; TktA 

Transketolase 2 isozyme, stationary phase induced/Non-oxidative Metabolism; TktB 

Translation initiation factor IF-2/Protein Translation; InfB 

Trigger factor/Cell division; TiG 

Tryptophanase^ryptophan Metabolism, Indole Production; TnaA 

Type II secretion pathway related protein^ransport; EtpE 

Tyrosine-protein kinase/Down regulates colanic acid production; WzC 

Uncharacterized protein/Hypothetical; ECs1547 

Uncharacterized protein/Hypothetical; ECs2891 

Uncharacterized protein/Hypothetical; ECs2991 

Putative oxidative stress defense protein/ Oxidative Stress; YggE 

Uncharacterized protein/Membrane Protein; YqjD 

UPF0042 nucleotide-binding protein/putative ATPase; YhbJ 

UPF0092 membrane protein, translocase^ransport; YajC 

UPF0337 protein, putative stress response protein/Osmotic Shock; YjbJ 

UPF0352 protein/Uncharacterized; YejL 

Uridylate kinase/Nucleotide Interconversion; PyrH 

Valine-tRNA ligase/Amino acyl tRNA Biosynthesis; ValS 



\Da, Kilodalton. 



o 
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SlyB, OmpA, BamA, BamD, TolC, OmpW, ElaB, YbjP, 
LppC, YqjD), chemotaxis (Trg), adherence (AidA-like, 
BssR, CsgG, CsgB, LpoB/YcfM, EspP), and cell division/ 
DNA replication (HrpB, ZapB, MukB, ParE, Ssb, Tig) 
(Table 4; Additional file 2: Table S2). 

Discussion 

This study provides a snapshot of various proteins 
expressed by 0157 in unfiltered, rumen fluid through a 
comparative analysis of the 0157 proteomic-profile in 
different media, growth conditions and incubation times. 
Interestingly, none of the reported (LEE, Shiga toxins) 
0157 virulence proteins were detected, under all con- 
ditions, in any media tested. Overall, fewer 0157 pro- 
teins were detected in more nutritionally complex 
RF-preparations versus LB and among these, differences 
were observed based on availability of oxygen, nutrients 
and incubation time. Also, the 0157-proteome in the 
RF-preparations included more proteins with diverse 
functions at 48 h than after 14 days of incubation. In 
fact, proteins associated with adherence, cell division and 
growth were identified only at 48 h. However, under all 
conditions, a selective expression of proteins with a role 
in cell structure, transport, metabolism, chemotaxis, mo- 
tility, resistance, stress and regulation was observed in 
RF-preparations , many of which were up-regulated in 
the unfiltered rumen fluid. The 0157 growth patterns 
and proteome expressed in the rumen fluid is suggestive 
of an adapting 0157, expending minimal energy, prepar- 
ing for survival and downstream intestinal colonization. 

Since adult cattle are often fed a maintenance diet 
with less protein until ready for feedlots, we decided to 
analyze 0157 growth dynamics in rumen fluid derived 
from animals on this diet. Rumen fluid from cattle fed a 
diet low in protein usually has a pH ranging from 6.2-6.8, 
and VFA concentrations at, 60-70% acetic acid, 15-20% 
propionic acid, 5-15% butyric acid [28-31]. The rumen 
fluid VFA and pH values were within the limits described 
for this diet for both animals used in this study (Tables 1 
and 2; 26-29). Irrespective of incubation times (14 days 
versus 48 h), 0157 exhibited very distinctive growth pat- 
terns in RF-preparations compared to LB. 0157 cultures 
in dRF, fRF and uRF were consistently at lower optical 
densities than LB, under both aerobic and anaerobic con- 
ditions. The anaerobic RF-preparation cultures never 
reached an OD 600 = 1.0 and the viable 0157 recovered 
were at substantially lower counts when compared to LB. 
The low OD readings and viable counts recovered from 
RF-preparation grown cultures may have been due to in- 
hibitory factors and /or limited nutrients in dRF, fRF, uRF, 
not seen in LB, having a bacteriostatic (aerobic) or bacteri- 
cidal (anaerobic) effect on 0157 and reflective of 0157 
growth in a stressful environment [11,32-36]. Using LB 
media for estimating viable counts may have helped recover 



the stressed bacteria [35]. Similar recovery of viable bacteria 
despite low OD reading has been reported among bacteria 
exposed to antimicrobial stress [36], and limited growth has 
been associated with bacteria entering into a stressed/ 
starved state or stationary phase [35-37]. 

Overall, fewer 0157 proteins were detected in RF- 
preparation cultures compared to LB, especially under 
anaerobic conditions. Irrespective of the media used to 
culture 0157, its anaerobic proteome was functionally 
associated only with cell structure, transport, metabol- 
ism, chemotaxis, motility, resistance, stress-related and 
regulation, and not 0157 virulence. Previous reports 
have demonstrated that 0157 virulence genes, especially 
the Shiga toxin and LEE-encoded genes, are down- 
regulated in LB compared to minimal media [38-40]. In 
addition, presence of trace amounts of glucose has also 
been shown to down-regulate LEE expression due to 
catabolite repression and/or acidic pH [38-40]. Hence, 
the lack of virulence gene expression in LB in this study 
conforms to those findings. Experiments with acid- 
stressed, starved bacteria have shown that these are 
likely to be more virulent only on recovery, and over 
time [35]. Even in minimal media that usually supports 
0157 virulence gene expression, several of these are 
suppressed as cultures reach the stationary phase [41]. 
Butyrate, a key environmental cue in LEE gene expres- 
sion was limited in the RF used in this study, which 
may have also caused the LEE suppression [9]. Condi- 
tioned media from unrelated cultures have been shown 
to suppress Shiga toxin gene expression while maintaining 
0157 growth or suppressing growth itself [33,35,42]. In 
fact, experimental studies have shown that it is easier to 
displace 0 157 in unfiltered rumen fluid versus autoclaved 
rumen fluid, by addition of "nonfermentable" sugars in the 
presence of the ruminal microflora [11]. Thus, the absence 
of 0157 virulence gene expression in RF-preparations 
may be reflective of the stressful growth environment, 
suppression due to nutrient limitations, lack of inducers, 
oxygen deprivation, pH fluctuations and inhibitory metab- 
olites released by resident microbiota. 

Previous studies have suggested development of acid 
resistance by Shiga-toxin producing E. coli (STEC) in 
the rumen as a means for better STEC survival through 
the stomach-like' acidic bovine abomasum [43,44] and 
have prescribed a role for glutamate-dependent acid 
resistance system (Gad system) and the tryptophanase 
(tnaA) enzyme toward this end [45]. Hughes et al, re- 
cently demonstrated that 0157 LEE expression is down- 
regulated while the Gad system is up-regulated in the 
rumen of cattle [46]. This observation made in animals 
being fed a grain diet, having a ruminal pH of 5.93, de- 
rived a role for the SdiA gene in sensing the acylhomo- 
serine lactone (AHL) signals in the rumen fluid and 
affecting differential expression of these genes. AHLs 
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formed by ruminal resident flora, are effective only 
under highly acidic pH and hydrolyze at neutral-alkaline 
pH [46,47]. Similarly, the Gad system that relies on the 
decarboxylation (gadA/B) of glutamate via proton con- 
sumption to increase cytoplasmic alkalinity is active at 
pH 4-4.6 [48]. However, other degradative amino acid 
decarboxylase and acid-resistance systems are activated 
in response to low pH (5.2 to 6.9), fermentative- 
anaerobic growth and stationary phase growth [48,49] 
and used more often than the Gad system to counter 
the deleterious effects of protons. We observed one such 
system, the arginine-dependent acid resistance system 
(Arg system) to be up-regulated in the RF-preparations 
after 48 h of anaerobic growth. Since the pH of the RF- 
preparations used in this study did not reach extreme 
acidic levels, the Gad system may not have been in- 
duced. In the Arg system, decarboxylation (speA) of 
arginine via proton consumption resulting in the forma- 
tion of agmatine stabilizes the cytoplasmic pH. Agmatine 
is either exported via the arginine- agmatine antiporter 
(aidC) or converted (speB) to putresceine as part of the 
polyamine biosynthetic pathway. 

Considering that 0157 is exposed to heat-shock, star- 
vation and stationary-phase-like growth in the rumen, it 
is possible that these factors enhance acid-tolerance in 
the bacteria through other mechanisms such as outer 
membrane changes and synthesis of proton transport- 
related protective proteins, as well [49,50]. Several stress 
(acid, low oxygen, osmolites, stationary phase) -responsive 
genes were expressed by 0 157 in this study, and included 
genes associated with the metabolism of arginine (speA, 
speB), lysine {lysU), formate (hyC), tryptophan (tnaA) and 
maltoporin (lamB), catalase (katG), DNA polymerase- 1 
(polA) and AidA-1 adhesin-like protein (aidA) [49-51]. 
Flagellar genes are differentially expressed under varying 
acid-stress conditions [51-53], and in our study, these 
genes were up-regulated in dRF and fRF but not uRF, sug- 
gesting less pH variation in the course of growth in uRF 
and limiting the role of flagella to motility alone. Stressed 
bacteria have been shown to be more adherent [35,40,53]; 
proteins associated with adherence (AidA-1 adhesin-like) 
and biofilm formation (BssR, CsgG, CsgB) were identified 
after 48 h incubation and not after longer incubation 
periods. Interestingly, several resistance' related proteins 
were up-regulated in RF-preparations, a subset of which 
(tellurite resistance, serine protease) have also been shown 
to contribute towards 0157 adherence [54,55]. This sug- 
gests that adherence may be critical during the initial 
phase of 0157 colonization and although LEE is sup- 
pressed, the bacteria rely on other mechanisms to adhere 
or form biofilms in the rumen. It has been observed that 
bacteria and protozoa in the rumen tend to adhere to the 
fibrous mat layers comprising of plant material to remain 
in the rumen and assist in the digestion of insoluble feed 



materials [56]. While this may not be in the case of 0157, 
initial adherence to or biofilm formation on available 
surfaces may give the bacteria time to adapt and survive 
the rumen environment [34]. It appears that much of 
the adaptive changes are initiated early in colonization 
as reflected in more stress-induced, structural integrity- 
related outer membrane proteins (AsmA, LptE, Lpp, 
NagA, SlyB, OmpA, BamA, BamD, TolC, OmpW, ElaB, 
YbjP, LppC, YqjD), and cell division and growth pro- 
teins, being expressed at 48 h. This supports the obser- 
vation that 0157 is maintaining slow growth in the 
RF-preparations as well. 

Conclusion 

Bottom-up proteomics provided a broad picture of dif- 
ferences in 0 157 protein expression after extended incu- 
bation in various media tested. Quantitative proteomics 
(iTRAQ) -based analysis of the 0157 anaerobic proteome 
expressed in uRF with all normal rumen flora was per- 
formed to more closely determine 0157 protein expres- 
sion in the bovine rumen. The cumulative results of all 
RF-preparation analysis suggested that rumen specific 
protein expression enables 0157 to adapt to this hostile 
environment and successfully transit to its colonization 
sites in the bovine GIT. To further verify our conclu- 
sions, we are evaluating the 0157 proteomic-profile as 
expressed in vivo in a rumen-fistulated cow, and con- 
firming the role of a subset of these adaptive' proteins 
in 0 157 survival. 

Additional files 



Additional file 1: Table SI. Bottom-up Proteomics Dataset. 
Additional file 2: Table S2. iTRAQ Proteomics Dataset. 
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